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We consider the Polyakov Nambu Jona Lasinio model with three massless quarks at high density 
and moderate temperature in the superconductive color flavor locking phase. We compute the 
critical temperature T c as a function of the baryonic chemical potential for the phase transition 
from the superconductive state to the normal phase. We find that T c is higher by a factor 1.5 -2 in 
comparison to the model containing no Polyakov loop. We also compute the specific heat C v near 
the second order phase transition and we show that the inclusion of the Polyakov loop does not 
change the value of the critical exponent. 



I. INTRODUCTION 



At small hadronic densities and sufficiently high temperature chiral symmetry is restored and the nature of the chiral 
phase transition can be investigated by various effective approaches. One of the most popular is the Nambu - Jona 
Lasinio (NJL) model [l[, describing the chiral transition in terms of the (qq) order parameter. At high temperatures 
one also expects a deconfinement transition Q. Its nature is rather clear in pure gauge theory, because, in absence of 
quarks, Quantum Chromodynamics at low temperature possesses a Z3 global symmetry, which is spontaneously broken 
at high temperature T. The order parameter for this phase transition is the Polyakov loop [3] whose expectation 
value vanishes in the disordered low temperature phase and is different from zero in the high T phase. 

The Polyakov loop is a SU(3) C matrix in color space given by (/3 = T _1 ) 

L(x) = Ve^fo dx 4 A 4 (x,x 4 ) _ h) 

For uniform A4 one gets an order parameter that can be written as follows in the Polyakov gauge 

<J> = _rp re i0(<t>3*3+<f>s*a) . (2) 
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for T — ► 00 one has $ = 1 and in the confined phase $ = 0. 

In presence of dynamical quarks there is no clear order parameter for the deconfinement transition because in this 
case the Z3 center of the SU(3) C gauge group is not a good symmetry. Though one cannot properly speak of a phase 
transition in this case, the T— dependence of the Polyakov loop can nonetheless be studied by numerical simulations 
and one still observes its rise from low to high temperatures on the lattice. 

An interesting and still debated 0] feature of these data is that chiral symmetry breaking and the decrease of the 
Polyakov loop occur at the same critical temperature It has been argued that a mixing between the effective 
models for chiral transition and the Polyakov loop dynamics might account for the approximate equality of these 
temperatures @ . A step forward has been obtained in Refs. Q and [1] , where the NJL model is studied in presence of 
a uniform extended gauge field A4. Its effect on dynamical quarks is obtained by identifying the parameters appearing 
in the Polyakov loop @ with an imaginary quark chemical potential. 

This modified NJL model (called Polyakov-NJL=PNJL) is characterized by a thermodynamical potential fl com- 
prising two terms, SInjl and U(T,<p). Qnjl contains the NJL thermodynamic potential modified by the inclusion 
of the imaginary quark chemical potential; in the mean field approximation it therefore depends on the chiral order 
parameter and on the Polyakov loop $. U(T,<fr) depends only on <E> and T and its parameters can be obtained by 
fitting pure gauge lattice QCD results. The PNJL model still stands on a conjectural basis. We do not analyze here 
its theoretical foundations. Nevertheless we do attempt to determine some of its possible physical implications. 

In Q the PNJL model has been extended to high baryonic densities for the case of two flavors (u, d) by including 
a quark chemical potential ji. At moderate fj, and small T a plausible model describing quark dynamics is the 2SC 
model [l(| dH characterized by condensation in the diquark antisymmetric color channel and decoupling of the strange 
quark. The two flavor approximation can only be valid at high, but not very high, densities. At these densities u and 
d quark masses can play a role, but their effect is included by considering also condensation in the qq channel. 

The aim of the present paper is to consider the case of higher densities, where all the three light quarks can form 
color superconductive pairs. The favored phase for sufficiently high density is the color flavor locking (CFL) state fl2| . 
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characterized by three massless quarks, qq condensation in spin 0, color and flavor antisymmetric state (for reviews see 
(l3|). This result was obtained in a NJL model, where the gluon interaction is mimicked by a four fermion interaction 
and one works in the mean field approximation. The dominance of the CFL phase can also be proved in QCD by 
way of one-gluon exchange; however this result is valid only at extreme densities (/i ~ 10 s GeV) [14]]. These densities 
are much larger than those presumably existing in the core of compact stars, where color superconductivity might 
be found. For these latter densities perturbative QCD is of little or no help. On the other hand the standard non 
perturbative method, i.e lattice QCD, is not applicable, as the quark determinant is complex at fi and MonteCarlo 
simulations are not directly usable. Therefore the four-fermion interaction remains as the only practical way to study 
the CFL phase. 

Though in this approximation gluon interactions are described by an effective four quark interactions, the Polyakov 
order parameter can nevertheless play a role, similarly to what happens for two flavors Q. The study of this role is 
the aim of this paper, where we present a preliminary study of the second order phase transition around the critical 
line T c {n). We will consider only the case of massless quarks, even though finite mass effects could be included either 
as free parameters or by considering condensation in both the diquark and quark-antiquark channels [HI ■ The reason 
for this neglect is that the inclusion of, say, the strange quark mass M s considerably complicates the analysis, because 
for M s ^ one should include electric and color chemical potentials to enforce electric and color neutrality. We will 
treat these effects, as well as the the extension to the gapless CFL (gCFL) phase [16j . in a future publication. 

Since we neglect mass effects, the free energy depends only on the order parameter $ and the unique gap parameter 
A (the role of the gap parameters due the symmetric color channels will be discussed below). Moreover we are 
interested only in the transition line between the CFL and the normal phase in the T — fx plane. Therefore we can 
use a Ginzburg-Landau (GL) expansion near the critical line. This approximation is discussed in Section [TTJ In 
Section Um we verify that the transition is continuous and compute the critical temperature as a function of the quark 
chemical potential fi. Our result is that the critical temperature T c is higher by a factor 1.5 - 2 in comparison with 
the treatment of CFL within the original NJL approximation. We also evaluate the critical exponent j3 that fixes the 
relationship between the gap parameter and the temperature near the phase transition and we find that including the 
Polyakov loop does not change the classical value (3 = 1/2. The discontinuity in the specific heat at the second order 
phase transition is also evaluated and a comparison of results obtained with and without Polyakov loop is performed. 
Finally, some concluding remarks are contained in Section [IVI 



II. THERMODYNAMICS OF THE THREE FLAVOR PNJL MODEL 



The model we study is described by the quark lagrangian 

£ = i> (i£>„y + /i 7o ) i> + C A ■ (3) 

In the above equation we have introduced the coupling of the quarks to a background temporal gauge field = 
gSftoA^Ta coupled to the quarks via the covariant derivative — — iA^] fi is the quark chemical potential. The 
term £a is responsible for color condensation. It can be obtained in the mean field approximation from a four fermion 
interaction term. In the CFL model one has 

A / \ SA 2 

£a = - y (vLts^^^'iC^. + h.c.) - — . (4) 

Eq. ((4]) describes the fact that in the ground state one has a non-vanishing expectation value of the di-quark field 
operator 

{^ai^ 1 ) OC ACapiCij! ^ . (5) 

The constant G in Eq. Q is the NJL four fermion coupling constant. In Eq. ([5]) we have neglected the color symmetric 
channel contribution, as one can prove that it becomes less and less important when one approaches the second-order 
phase transition. 

Once the quark lagrangian is specified, the mean field free energy of the CFL quark matter is easily obtained by 
integration over the fermion fields in the generating functional, namely 

V=U{T,4>) + — -2 Tl 22J (2^ l0 H f J ' (6) 

where to n — irT(2n + 1) are the fermion Matsubara frequencies, and is the inverse fermion propagator in the 
mean field approximation, whose explicit form can be found in Ref. [16] . S* -1 is in principle a 72 x 72 matrix in color, 
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flavor, spin and Nambu-Gorkov indices. In the high density limit the effect of the antiparticles can be neglected; 
moreover, one can split the left-handed and the right-handed quark contributions to the free energy, since the quarks 
are massless and the condensation does not mix quarks with opposite chirality. Thus S^ 1 is reduced to a 18 x 18 
matrix. It can be rearranged to a block diagonal form, with a 6 x 6 matrix describing the propagation of u r ,d g ,Sb 
quarks, and three 4x4 matrices describing the propagation of d r ,u g , and s r ,Ub, and s 9 ,(4 quarks. This allows a 
straightforward extraction of the quasiparticle dispersion laws, much in the same way as in the analogous evaluation 
contained in [15], the difference being that there Ai = 0, A2 = A3 and here Ai = A2 = A3 = A. 

In Eq. ^ we have introduced the part of the thermodynamic potential IA(T, (ft) which describes the dynamics of 
the Polyakov loops in absence of dynamical quarks. In principle various forms can be used 0, Q; for definiteness we 
adopt the form proposed in @ 

U(T, <fi) = T 4 { - + b(T) ln[l - 6$*$ + 4($ 3 + $* 3 ) - 3($*$) 2 ] } (7) 



where 

a(T) = a + 01 ) + «2 {jr ) , b ( T ) = b * [f ) ' ( 8 ) 

The Polyakov loop $ can be expressed in terms of one parameter <j> = ^3, as the other parameter 0g can be always 
be absorbed by a redefinition of ^3 . It is given by 

1 + 2 cos( 8(h) , , 

$ = $* = — — KHVJ . (9) 

Numerical values of the coefficients have been fitted in using lattice data [l8| : 

a = 3.51, ai = -2.47, a 2 = 15.2, 63 = -1.75 , (10) 

together with the dcconfinemcnt temperature To = 270 MeV . The use of a definite form for U(T, </>) is not a limit 
of our computation because other functional dependences produce similar results as they are derived from the same 
lattice data set. The NJL coupling G should in principle depend on $ because the four fermion coupling is induced 
by gluon dynamics. However following [6( we will neglect this effect. In final results we will trade G for the value Ao 
of the CFL gap at T = using the weak coupling formula [l!| 



3 _ 6a* 2 l 2(5 
G 7T 2 2^A 



2 ln TT— ■ (11) 



In the above relation S is an ultraviolet cutoff, introduced to ensure ultraviolet convergence of the loop integrals. By 
means of Eq. (jTTJ) the coupling strength is represented by the parameter Ao. 
Performing the summation over the Matsubara frequencies in Eq. ([5]) one gets 

fl = U(T^) + — -2£ J j^{T\n 1 + + - p + M ) + 

+ rin|l + e-^| +»(£,• -p-fi)) ■ (12) 



In Eq. (fT2"j) €j are the energies of quasiparticles and ej are obtained from €j by the substitution £ = p — \i —> p + fi. The 
terms with ij correspond to antiparticles. They are put here for completeness but omitted in numerical evaluations. 
The dispersion laws for the nine quasiparticles can be derived by the standard methods. Since the resulting expressions 
are cumbersome and our study is limited to the critical line we present here their expression only for small values of 
the gap parameters, ie in the GL approximation. One obtains 

, / , A 2 8f 2 + 2 . A 2 

ei = e 2 = K + W 



= a 1 



2£ 4f 2 + 4> 2 4£ 2 

4A 2 



4^ 2 + (j) 2 
A 2 

e 4 = 4 = C + i^+ — 
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e 6 = £s = £ + i<t 



A 2 



(13) 



We have also computed the coefficients of the C(A 4 ) term but we do not report them here (they are needed to control 
that the phase transition is continuous at T c and to compute the gap, see below). 

The gap parameter A and the background gauge field at a fixed temperature and chemical potential are obtained 
solving the equations 



0, 
0. 



dA 

Near the critical temperature T c one can expand f2 in Eq. (fT2|) as follows 

fi(A,0)-f7(O,0)^|A 2 + ^A 4 . 



(14) 
(15) 

(16) 



The critical temperature T c is obtained at a fixed /i as in the usual BCS theory by solving the equation a(T c ) = 0, 
with a given by 



6 

a = G 



ct 2 ^ v f V [3 ■ (p - e? + 4> 2 ■ {ii + 3a] 

3 ^ V J -5 (*o - [Co + 2 + 2 ] l(lo - 2 + </> 2 ] 



lo—iu n 



In- 



25 
23A 



1 

3 ./_ 



(17) 



and / defined by 



2e 



4C 2 



— tanh — 2^——— 

6 2 2 4£(2£ - i<j 



■ tanh 



m - i<t>) 



(18) 



This expression for a is identical to the result obtained by (|12p using the dispersion laws (fT3"|) up to A 2 . On the other 
hand the the coefficient j3 is given by 



= T 



2tt 2 



/■ 5 



8A 4 -^(/ ,e,$) 



^ {(ig-£ 2 )[(io+0 2 + Wo-0 2 + W 



iQ—iuj-n 



(19) 



where 



^(/ , £, *) = 62 8 - Z 6 • (24£ 2 + 5^ 2 ) + l 4 Q ■ (36£ 4 + 21^0 2 - 40 4 ) 

-Z 2 , • (24f + 26£ 4 2 - 24£ 2 4 + 6 ) + £ 2 • (6£ 2 - </> 2 ) • (£ 2 + 2 ) 2 



(20) 



The summation over Matsubara frequencies in the expression of (3 can be performed analytically, but the final 
expression is involved and we omit it for simplicity. 



III. NUMERICAL RESULTS 



To get the critical temperature T c we solve the equation a(T c ) = with a given by Eq. (|17[) and <f> obtained by Eq. 
(pjj. We have checked that the phase transition is of the second order since (3(T C ) > 0. It is well known that in the 
case 4> — one has in the CFL phase T c /Ao ~ 0.71, see for example [2fj- On the other hand in the 2SC phase one 
has T c /Aq ~ 0.57 as in ordinary BCS superconductors. This difference is related to the fact that in the CFL model 
one has eight gapped modes with gap A and one mode with gap 2 A. 

The result of the numerical evaluation of T c is shown in Fig. [1] where we plot the ratio T c /Aq with (solid line) and 
without (dashed line) Polyakov loop, at the reference value fi = 500 MeV. We notice that introducing self-consistently 
the parameter <f> implies a significant increase of the critical temperature. This effect has been noticed also in the two 
flavor model 0]. 
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FIG. 1: Ratio T c /Ao against Ao (MeV), with (solid line) and without (dashed line) Polyakov loop, at /i = 500 MeV. 



Next we turn to the behavior of the gap parameter A for temperature close to T c . We find 



A(T) 



/•■(A„)(l- J 







(21) 



with (3 — 1/2. The value of the critical exponent is the same as in BCS superconductors. However the presence of 
the Polyakov loop affects the constant k in two ways. First, it gives it a dependence on Ao that is absent in the BCS 
and in the two flavor color superconductor. Second, it changes its numerical values. For example for the 2SC case 
k ~ 3.1; in the present case k — 1.7 and 2.2 for A = 40 and 100 MeV respectively. 
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FIG. 2: Specific heat C v (Units: 10 7 MeV 3 ) against T/T c , with (solid line) and without (dashed line) the Polyakov loop, for 
H = 500 MeV and A = 25 MeV. 

The knowledge of A(T) near T c allows to determine some thermal properties of the model. For example we compute 
the specific heat as a function of the temperature, near T c . It is given by 

C. = -t£^*> . (22, 

We show the result of this calculation in Fig. [21 with (solid line) and without (dashed line) Polyakov loop, for 
fi = 500 MeV and Ao = 25 MeV (for other values of Ao we find qualitatively similar results). We notice that 
including the Polyakov loop slightly decreases the specific heat and increases a bit its discontinuity around T c . 
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IV. CONCLUSIONS 



In this paper we have studied the effect of the inclusion of the Polyakov loop on the NJL description of the CFL 
model. We have restricted our attention to a temperature range close to the critical temperature of the second order 
phase transition. We have found that introducing the Polyakov loop significantly increases the critical temperature, 
the effect being more important in the weak coupling regime. This increase may have some phenomenological conse- 
quences, both for astrophysical systems and for future experiments at GSI, if the proposed facility SIS100/200 [2l[ 
will be able to reach the hadronic densities needed for color superconductivity. Needless to say, one has to stress the 
heuristic use of the Polyakov loop when quarks are dynamical. Already their presence destroys the center symmetry 
of pure gauge QCD. More theoretical investigation will be needed on the PNJL model to ascertain its possible regions 
of validity. Nevertheless we felt that it is useful to investigate the effect of the Polyakov loop in some portions of the 
QCD phase diagrams where a direct QCD treatment is not available at the present. 

We have studied the behavior of the gap parameter A(T) for T w T c , showing that the Polyakov loop does not 
modify the critical exponent [3 — 1/2, but only the pre-factor. In ordinary superconductor the pre-factor does not 
depend on A ; on the other hand, the presence of the Polyakov loop results in a pre-factor dependent on the strength 
of the coupling. 

A quantity of interest is the specific heat C v since it can be measured experimentally. At the second order phase 
transition C v is discontinuous, in the superconductive phase being larger than in the normal phase. Although the 
effect of the Polyakov loop is to decrease the absolute value of C v , the discontinuity AC V with <fr ^ is larger than 
the corresponding value at <p = 0. 

Further developments include the treatment of the strange quark mass, as well as the study of the thermodynamics 
of the CFL superconductor with Polyakov loop at lower temperatures. While the effect of the Polyakov loop is not 
expected to modify strongly the thermodynamics at small T, it is known that the finite strange quark mass affects 
the phase diagram of QCD and its role in the present model may be of interest as well. 
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